Jet evolution in Yang-Mills- Wong simulations 



Bjorn Schenke' 



"Department of Physics, McGill University, Montreal, Quebec, H3A 2T8, Canada 



Abstract 



We present results for collisional energy loss and momentum broadening of high momentum 
partons in a hot and dense non- Abelian plasma obtained by solving the coupled system of Yang- 
Mills-Wong equations on a lattice in real time. Including hard elastic collisions among the parti- 
cles we obtain cutoff independent results for the collisional energy loss dEjdx and the transport 
coefficient q. The latter is found to receive a sizable contribution from a power-law tail in the 
transverse momentum distribution of high-momentum partons. We further argue that the effect 
of instabilities on jet broadening should be accessible by experiment when employing jet cones 
with elliptical bases or studying correlations within the cone in full jet reconstruction methods. 



1. Introduction 

High transverse momentum jets produced in heavy-ion collisions represent a valuable tool 
for studies of the properties of the hot parton plasma produced in the central rapidity region HI. 
We study the collisionsal energy loss and momentum broadening of high momentum partons 
employing a numerical simulation of the Boltzmann-Vlasov equation, which is coupled to the 
Yang-Mills equation for the soft gluon degrees of freedom. Soft momentum exchanges between 
particles are mediated by the fields, while hard momentum exchanges are described by a collision 
term including binary elastic collisions. This way, we are able to provide an estimate of the 
coupling of jets to a hot plasma which is independent of infrared cutoffs. 

We also simulate plasmas with local momentum-space anisotropies, which occur in heavy- 
ion collisions due to the longitudinal expansion during the very early stages of the plasma evo- 
lution. These plasmas develop Chromo-Weibel instabilities [2], which lead to the formation of 
long-wavelength chromo-fields with E z > B z and B ± > E ± . The strong fields affect the propaga- 
tion of hard partons, leading to an asymmetry of the jet shape in rapidity Arj and azimuthal angle 
A<p, which should be accessible using new jet reconstruction measurements |0]. 

2. Boltzmann-Vlasov equation for non-Abelian gauge theories 

We solve the classical transport equation for hard gluons with SU(2) color charge including 
hard binary collisions 



where / = f(x,p,q) denotes the single-particle phase space distribution. It is coupled self- 
consistently to the Yang-Mills equation for the soft gluon fields 



P» [d M + gq a F; v d v p + gf^A^x)^] f = C, 



(1) 




(2) 
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with = (l,p/p). The collision term contains all binary collisions, descibed by the leading- 
order gg — > gg tree-level diagrams. 

We replace the distribution f(x,p,q) by a large number of test particles, which leads to 
Wong's equations J3] 

x,(f) = V,(?) , 

Pi(t) = ?9 ;(0(F(0 + v ( (0xB"(r)), 

qM = -ig^(t)[A,(t\ qi (t)], (3) 

for the z'-th test particle, whose coordinates are x,(f), P;(0> and q"(t), the color charge. The time 
evolution of the Yang-Mills field is determined by the standard Hamiltonian method |5[] in A = 
gauge. See Ual7llMl2lllO(] for more details. The collision term is incorporated using the stochastic 
method Hill . The total cross section is given by <T2-^2 = f k ll 2 2^dq 2 , where we have introduced 
a lower cutoff k*. To avoid double-counting, this cutoff should be on the order of the hardest 
field mode that can be represented on the given lattice, k* ^ n/a, with the lattice spacing a. This 
way soft momentum exchanges are mediated by the fields while hard momentum exchanges are 
described by the collision term. For a more detailed discussion see ifioll . 



3. Jet broadening in an isotropic plasma 

We first consider heat-baths of particles with different densities n g and temperatures T. We 
simulate an undersaturated gluon plasma due to numerical restrictions (see tlOfl ) and later ex- 
trapolate to thermal systems with physical densities and temperatures. For a given lattice (or k*) 
we take the initial energy density of the thermalized fields to be J d 3 k/(2n) i kfB 0se (k)®(k* - k), 
where fsoseik) = n g /(2T 3 ({3))/(e k ^ T - 1) is a Bose distribution normalized to the assumed parti- 
cle density n g , and ( is the Riemann zeta function. The initial field spectrum is fixed to Coulomb 
gauge and A,- ~ l/k. 

We study colorless bunches of test particles that allow us to restrict to collisional energy loss 
and momentum broadening due to elastic collisions only. Both 

and the differential energy loss dE/dx are independent of the separation scale k* . A is the mean 
free path and p ± denotes the momentum transverse to the initial jet momentum. Fig. Q] depicts 
the contributions to q and dE/dx due to soft and hard collisions separately, as well as the total. 

Studying the dependencies of q and dE/dx on the density, temperature, and hard parton 
energy reproduces results from perturbative QCD (see [10] for a detailed analysis). Using these 
dependencies we can extrapolate q and dE/dx to temperatures around 400 MeV. Adjusting the 
color factors as appropriate for SU(3) and extrapolating to the thermal density of gluons we find 
q ~ S.e + OJGeV^nT 1 and dE/dx * 1.6 + 0.4 GeV for 1 (E = 19.2 GeV). 

Next, we present the full /^-distribution of the high-momentum partons traversing the hot 
medium in order to assess the relative contributions from various processes to its first moment 
q. We find that over time the initial ^-function broadens to a Gaussian distribution with a power- 
law tail. Fig. |2] shows the distribution of the high-momentum test particles after t a 5.2 fm in 
a double-logarithmic plot versus p]_/t. The low-p ± part follows a Gaussian distribution in p ± . 
The power-law tail at large p ± behaves approximately as p~ 4 . This is expected for particles 
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Figure 2: p\ -distribution of the high-momentum (E/T = 
48) partons after t a 5.2 fm for T = 4GeV and n = 
20 fnT 3 (q x 5.16GeV 2 fnr'). 
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Figure 3: Nuclear modification factor Raa(P±) of jets due 
to elastic energy loss in a classical Yang-Mills field pro- 
duced in the early stage of a relativistic heavy-ion colli- 
sion at RHIC. The band indicates the uncertainty due to 
the extrapolation of dE/dx to physical temperatures. 



experiencing only few scatterings since in the high-energy limit the differential cross section 
dcrldp\ ~ p'] 4 . We also indicate the value of q to show that the power-law tail contributes 
significantly to this transport coefficient. 

We also provide an estimate for the nuclear modification factor Raa due to elastic energy 
loss in a classical Yang-Mills field. This is of relevance for collisions of heavy nuclei at high 
energies: the large number of gluons produced in the central rapidity region can be described as 
a classical field for a short time (see e.g. HH) until the field modes decohere and thermalize @]. 
From the simulations presented above, we obtained the elastic energy loss dE/dx. For large k* 
(on the order of the "saturation momentum" Q s ) most of the energy density is due to the classical 
field. In a thermal system at T — 400 MeV the energy density is about 17 GeV/fm 3 , which is 
an appropriate average over the first 1 fm/c of a central Au+Au collision at RHIC energy [12]. 
An estimate for the nuclear modification factor Raa at the parton level (neglecting hadronization) 
can be written as lfl3l[loll 

dNf/d 2 p±dy 

Raa{ P^ = J , J T = (1 " «PJf ■ (5) 
dNj/d / -p±dy 

3 



e denotes the fractional energy loss up to t — 1 fm and dNj/d 2 p±dy ~ l/p" +2 is the initial p ± 
distribution of jets, where n w 4. Fig. [3]shows the result for Raa- The experimentally observed fiat 
Raa ~ 0.2 can not be fully accounted for by early-stage elastic energy loss in the classical field 
background, but our result shows that this contribution is significant and can not be neglected. 



4. Jet broadening in an unstable plasma 

As studied in detail in J9(], instability growth in unstable plasmas leads to a direction depen- 
dent q. We describe the broadening of the hard (now colored) test particles in the transverse ± 
and longitudinal z directions via the variances q ± := ^((A/? ± ) 2 ) ,q z := ■j l {(Ap z ) 2 ) . The ratio 
q z /q± can be roughly associated with the ratio of jet correlation widths in azimuth and rapidity: 
yk~zlWx ~ (Ari)/{A<f>). The numerical simulations of unstable plasmas obtain 19[] 

<Ai/>/<A0>*L5. (6) 

We propose to look for this effect by using the greatly improved experimental methods for full 
jet reconstruction J5, 14, 15], In particular, one should use cones with an elliptical base and 



determine its major axis in the cf>-rj plane by maximizing the energy contained in the cone. If 
the system is isotropic, the angle should be uniformly distributed while the effect of instabilities 
should lead to a preferred orientation of the major axis in the rj direction. Another, possibly more 
practical way of experimentally determining such an anisotropy of the jet is to keep the circular 
base and study correlations among the particles within the cone. 
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